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Abstract 

Hierarchy of one-parameter families of chaotic maps with an invari- 
ant measure have been introduced, where their appropriate coupUng 
has lead to the generation of some coupled chaotic maps with an in- 
variant measure. It is shown that these chaotic maps (also the coupled 
maps) do not undergo any period doubling or period-n-tupling cascade 
bifurcation to chaos, but they have either single fixed point attractor 
at certain values of the parameters or they are ergodic in the com- 
plementary region. Using the invariant measure or Sinai-Rulle-Bowen 
measure the Kolmogrov-Sinai entropy of the chaotic maps (coupled 
maps) have been calculated analytically, where the numerical simula- 
tions support the results. 

Keywords: coupled chaotic map, chaos, Invariant measure, 
Kolmogrov-sinai entropy, Lyapunov characteristic exponent, 
Ergodic dynamical systems. 
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In recent years chaos or more properly dynamical systems have become an 
important area of research activity. One of the landmarks in it was the in- 
troduction of the concept of Sinai-Ruelle-Bowen (SRB) measure or natural 
invariant measure. This is roughly speaking a measure that is supported on 
an attractor and it describes the statistics of the long time behavior of the or- 
bits for almost every initial condition in the corresponding basin of attractor. 
This measure can be obtained by computing the fixed points of the so called 
Frobenius-Perron(FP) operator which can be viewed as a differential-integral 
operator, hence, exact determination of invariant measure of dynamical sys- 
tems is rather a nontrivial task, such that invariant measure of few dynamical 
systems such as one-parameter families of one- dimensional piece wise linear 
maps like the Baker and tent maps [0, 0, ^, Q or unimodal maps like logistic 
map for certain values of its parameter ( Ulam-Von Neumann map), one di- 
mensional cubic and quintic maps ^ can be derived analytically. In most 
of cases only numerical algorithms, as an example Ulam's method [0, 
are used for computation of fixed points of FP operator. 
Here in this article, we give a hierarchy of one-parameter families $Ar(Q;, x) of 
maps of interval [0, 1] with an invariant measure, such that Ulam-Von Neu- 
mann map, one dimensional cubic and quintic maps correspond to three first 
maps of this hierarchy for particular value of their parameter. Also the map 
appearing in the study of field driven thermostated systems [Q is topolog- 
ically conjugate to the first map ( generalized Ulam-Von Neumann) of our 
hierarchy of family of one parameter chaotic maps. Then, by an appropriate 
coupling of these maps, we generate some coupled chaotic maps with an in- 
variant measure. We have been able to derived analytically their invariant 
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measure. Using this measure, we have calculated analytically, Kolmogorov- 
Sinai(KS) entropy or equivalently positive Lyapunov characteristic exponent 
of these maps, where the numerical simulations support the analytic calcu- 
lations. Also it is shown that these maps have another interesting property, 
that is, they bifurcate to chaotic regime without having any period doubling 
or period-n-tupling scenario. 

The paper organized as follows. In section II the hierarchy of one parame- 
ter families of chaotic maps have been introduced. Section III is devoted to 
the appropriate coupling of these maps. In section IV we talk about their 
invariant measure. In section V, using the invariant measure, we give their 
KS entropy. Finally in section VI we try to support the analytic calculations 
with the numerical simulations. The paper ends with brief conclusion and 
five appendices which contains all algebraic calculations and proofs. 

II Hierarchy of one- parameter families of chaotic 
maps 

One of the exactly solvable one dimensional maps or the map with invariant 
measure is the well known Ulam-Von Neumann map 

f{x) = x{l — x) (Logistic map for /i = 1). In reference ^ the Ulam-Von 
Neumann map with KS-entropy ln(2) has been generalized to exactly solvable 
one dimensional cubic and quintic maps with KS-entropy ln(3) and ln(5), 
respectively. Here we generalize these maps to a Hierarchy of one parameter 
families of maps with invariant measure, where the map of Reference is 
topologically conjugate to the first map (generalized Ulam-Von Neumann) of 
this hierarchy. 
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These one-parameter families of chaotic maps of the interval [0, 1] with an 
invariant measure can be defined as the ratio of polynomials of degree N: 



.(1), a2(^l + (-l)^Fi(-iV,iV,i,x 



(«2 + 1) + («2 _ l)(-l)N^F,i-N,N, l,x) 



2' 

,.2(rr ( /-:\\2 



fII-11 



l + (a2-i)(T;v(v^)2) ' 
where is an integer greater than one. Also 

2F^{-N,N,^,x) = (-l)^cos(2Ararccosv^) = (-l)^T2,v(v^) 

is the hypergeometric polynomials of degree N and T/v are Chebyshev poly- 
nomials of type respectively. Obviously these map the unit interval 
[0, 1] into itself. 

^N{a,x) is (N-l)-model map, that is it has {N ~ 1) critical points in unit in- 
terval [0, l](see Figures 1,2), since its derivative is proportional to derivative 
of hypergeometric polynomial 2Fi{—N,N,^,x) which is itself a hypergeo- 
metric polynomial of degree (A^ — 1), hence it has (A^ — 1) real roots in unit 
interval [0, 1]. Defining Shwarzian derivative S{^i\f{x,a)) as[|l^]: 

^'j^{x,a) 3f^'l^{x,a)y f ^'l^{x,a)\' 1 f ^%{x,a)Y 



<l>'jv(a;,a) 2\^'j^{x,a)J \^'j^{x,a)J 2\^'j^{x,a) 



with a prime denoting a single differential, one can show that: 
S {^n{x, a)) = S (^2Fi{-N, N, \.x)'^ < 0, 



since -^{2Fi{—N,N, can be written as: 



Xi 



with < a;i < a;2 < ^3 < .... < xn-i < 1, then we have: 



1 \ — 1 1 /N-i ^ 2 
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Therefore, the maps $iv(x, a) have at most + 1 attracting periodic orbits 
]T2| . As we will show at the end of this section, these maps have only a single 



period one stable fixed points or they are ergodic. 

Denoting n-composition of functions ^(x, a) by it is straightforward to 
show that the derivative of at its possible n periodic points of an n-cycle: 
X2 = ^(xi, a), X3 = $(x2, a), ■ ■ ■ , xi = a) is 

n 

I 1=1 A($(i>2)o$(M)o^.o$(i'2)(x,a)yi= n I 

(XiOC diCC I ^ 

fc=l 

(II-2) 

since for G [0, 1] we have: 

mini^a^ + (1 — a^Xk)) = min{l, a^), 

therefore, 

mm — = — rmnil^a ) 
dx \a J 

Hence the above expression is definitely greater than one for < a < A^, that 
is: maps do not have any kind of n-cycle or periodic orbits for < a < A^, 
actually they are ergodic for this interval of parameter. From (II-2) it follows 
the values of | \ at n periodic points of the n-cycle belongs to interval 

[0,1], varies between (Na)"' and (■^)" for a < and between (^)" and 
(Na)"' for a > N, respectively. 

From the definition of these maps, we see that for odd N, both x = and 
X = 1 belong to one of the n-cycles, while for even N, only x = 1 belongs to 
one of the n-cycles of $Ar(x, a). 

For a < (-^) {a > N), the formula (// — 2) implies that for those cases in 
which X = (x = 1) belongs to one of n-cycles we will have | -^^^^^ |< 1, 
hence the curve of starts at x = (x = 1) beneath the bisector and 
then crosses it at the next (previous) periodic point with slope greater than 
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one, since the formula {II — 2) implies that the slope of fixed points increases 
with the increasing (decreasing) of | |, therefore at all periodic points of 
n-cycles except ior x — {x — 1) the slope is greater than one that is they 
are unstable, this is possible only if x = (x = 1) is the only period one 
fixed point of these maps. 

Hence all n-cycles except for possible period one fixed points x — and x — 1 
are unstable, where for a e [0, j^] , the fixed point x = is stable in maps 
^j\f{x,a) (for odd integer values of N), while for a G [N, oo) the a; = 1 is 
stable fixed point in maps (see Figures 3,4 for N — 2,3). 
Below we also introduce their conjugate or isomorphic maps which will be 
very useful in derivation of their invariant measure and calculation of their 
KS-entropy in the next sections. Conjugacy means that the invertible map 
h{x) — maps / = [0,1] into [0, oo) and transform maps $jv(x, a) into 
$jv(2^, ce) defined as: 

^n{x, a) ^ ho ^n{x, Oi) o h^~^^ = tan^(A^ arctan V^) (II-3) 

As an example we give below some of these maps: 

a\2x-lf 
4a;(l - x) + a^(2x - 1) 



^2(«, x) = r. JL ..2 ^ (n-4) 



a'^x{Ax — 3)^ 
" a2^(4x-3)2 + (l-x)(4x-l)2' ^ 

, , - Sx{l - x)f 

<P4i«, X) ^2(1 _ 8^(1 _ ^))2 + 16^(1 -x)(l- 2xy ' ^ ^ ^ 

^ / X q;2x(16x2 -20x + 5)' 

$5(q;,x) = ^ 5 5 . (II-7) 

^ ' a^x (16x2 - 20x + 5)^ -X (16x2 -20x + 5)^1 ^ ' 



III Hierarchy of coupled map 

Using the hierarchy of families of one-parameter chaotic maps (//—I), we can 
generate new hierarchy of coupled maps with an invariant measure. Hence, 
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by introducing a new parameter g as a coupling parameter we form coupling 
among the above mentioned maps, where they can be coupled through (3 and 
a functions defined as: 

(3{x) = {J^o + gxf, 



(3{X) 



X) 



(III-l) 
(III-2) 



with Bn{x) and Ai>f{x) defined as: 



k=0 



k=0 

Now, the hierarchy of the coupled maps can be defined as: 

^1 ^Ni{Xl,ai{x2)) 

Their conjugate or isomorphic maps are defined as: 



(III-3) 



(III-4) 



Ni,N2 



— < 



(III-5) 



Xn^ {xi, ai{x2)) ^ho Xn^ oh^ ^) = tan2(A^i arctan y/x{) 

Xn2{x2, 0C2) — ho Xn^ o — 4j tan^(A^2 arctan ^/x^). 

(III-6) 

Below we give some of coupled map as an illustration 



^2,2 = < 



^1 



aj(x2)(2xi-l)^ 
4xi(l-xi)+a^(x2)(2x'i-l)2 



«i(2a;2-l)' 



(III-7) 



I - 4x2(l-a;2)+a^(2x2-l)2 



and 



^1 



2,3 



Q:f(x2)(2xi-l)^ 
4xi{l-xi)+ai{x2){2xi-l)^ 



(III-8) 



Y '-'2-' 

. ^2 — a|a;2(4a;2-3)2+(l-i;2)(4x2-l)2 
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with 



ai{x2) = 



l+f3{x2) \ (3{X2)' 



P{X2) 



ai 



+ gxf 



V 2 - ai 

One can show the map $2,3, is ergodic, for < ai < 2 and ^ < ^2 < 3, while 
for 2 < q;i < 00 and < q;2 < | (3 < q;2 < 00) it has fixed point at [0, 0] 
([0,1]) (see Fig. 5). 

a'j(x2)xi(4:Xi-3)^ 
;a;2)xi(4a;i-3)2+(l-a;i)(4a;i-l)2 

$3 9 = < (III-9) 



X 



2 4x2(l-a;2)+a^(2a;2-l)2 



Similarly, one can show that the map $3^2, is ergodic for | < ai < 3 and 
< q;2 < 2, while for < ai < | (3 < ai < 00) and < q;2 < 2 
(2 < q;2 < 00) it has fixed point at [0, 0] ([1, 0]) (see Fig. 6). 

^1 ~ a2(a;2)o;i(4a:i-3)2+(l-o;i)(4o;i-l)2 
$3 3 = < (III- 10) 



a:|a;2(4o;2— 3)2 



aSo:2(4o:2-3)2+(l-a:2)(4o:2-l)2 



with 



ai{x2) 



(3{x) 



3/^(0:2) + 1 /3(X2) 
3 + /3(x2) ^ I3{X2)' 

gxf. 



V 3q;i — 1 

Also the map $3,3, is ergodic for | < cti < 3 and | < q;2 < 3, while for 
< CKi < I (3 < CKi < 00 ) and < q;2 < | (3 < q;2 < 00 )it has fixed point 
at [0,0] ([1,1]) (see Fig. 7). 



^1 



a^(o:i)(l-8o:i(l-3:i))2 

a^(3;i)(l-8o;i(l-3;i))2+162;i(l-o:i)(l-2x2)2 



4,3 



(III-ll) 



a^a:2(4a:2-3)2 

a2a;2(4o:2-3)2+(l-o:2)(4o:2-l)2 
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and 



^4,4 = < 



alixi)il-8xi{l-xi)f 



a^(xi){l-8a:i(l-xi))2+16xi(l-xi)(l-2z2)2 



a2(l-8a;2(l-a;2))2 



with 



. 2 Q2(l-8x-2(l-a;2))2+16a;2(l-a:2){l-2x2)2 
/?2(X2) + 6f3{x2) + 1 



ai(x2) 



PiX2) 

/9(x2)(4/3(x2)+4) ^ /3(X2)^ 



4 - 6ai - 4./1 - 2ai + 2a{ 



(III-12) 



2(«i-4) 

The maps $2,2, '^'4,3 and $4^4 are also either ergodic or they only single fixed 
points at [0, 0] or [0, 1]. 



IV Invariant measure 

Dynamical systems, even apparently simple dynamical systems which are de- 
scribed by maps of an interval can display a rich variety of different asymp- 
totic behavior. On measure theoretical level these type of behavior are de- 
scribed by SRB [ p!3| , |T^ or invariant measure describing statistically station- 



ary states of the system. The probability measure /i on [0, 1] is called an 
SRB or invariant measure of the map y = ^n{x,C() given in (// — 1), if it 
is $Ar(a;, a)-invariant and absolutely continuous with respect to Lebesgue 
measure. For deterministic system such as <l>Ar(x, a;)-map, the $Ar(x,a)- 
invariance means that, its invariant measure fulfills the following formal 
FP integral equation 

This is equivalent to: 
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defining the action of standard FP operator for tlie map $Ar(a;, a) over a 
function as: 

We see that, the invariant measure //(x) is actually the eigenstate of the FP 
operator corresponding to largest eigenvalue 1. 

As it is proved in Appendix A, the invariant measure iJ>^j^{x,a){x, has the 
following form: 

pi{x) = - , ^ (IV-3) 

with /3 > 0, is the invariant measure of the maps ^n{x, a) provided that, we 
choose the parameter a in the following form : 

r (iv-l) i 

a = ^^=0 , (IV-4) 

where the symbol [ ] means greatest integer part. 

As we see, the above measure is defined only for > 0, hence from the 
relations {IV — 3), it follows that the maps $Ar(x, a) have invariant measure 
only for a G (-^, N) for odd value of N, and for even value a G (iV, oo). For 
other values of a these maps have single attractive fixed points, which is the 
same as the prediction of the previous section. 

Similarly the probability measure ii for coupled chaotic maps ^Ni,N2 given 
in (/// — 4) fulfills the following formal FP integral equation 

A*(?/i,l/2) = / dxi j dx2S{yi- X'^{xi,X2))d{y2- X2{xi,X2))n{xi,X2), 
Jo Jo 

which is equivalent to: 

Kyi,y2) ^ E \ J{Xl,X2) \l^{Xi,X2), (IV-5) 
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where J{xi,X2) is the Jacobian of the transformation ^n-^^n^ which is equal 
to 

again defining the action of standard FP operator for the map $Afi,Ar2 o'^^r a 
function as: 

P'^.,,.Jiy)= E fix)Jix^,X2). (IV-6) 

We see that, the invariant measure yu(xi, X2) is also the eigenstate of the FP 
operator -P$jv^ corresponding to largest eigenvalue 1. 
As it is proved in Appendix B the invariant measure ^^(xi,X2) has the 
following form: 

1 y% ^ 1 VA(a^2) 

^ v/2;2(l-X2)(/52 + (1 - P2)x2) ^ ^xi(l-a;i)(A(x2) + (1 - /3i(a;2))a;i) ' 

(IV-7) 

with /?2 > and (3i{x) > given in {III - 1). 



V Kolmogrov- Sinai entropy 



KS-entropy or metric entropy measures how chaotic a dynamical system 
is and it is proportional to the rate at which information about the state 
of dynamical system is lost in the course of time or iteration. Therefore, 
it can also be defined as the average rate of information loss for a discrete 
measurable dynamical system ($7v(a^, «), /w)- By introducing a partition a = 

Ac{ni, n^) of the interval [0, 1] into individual laps Ai one can define the 

usual entropy associated with the partition by: 

71(7) 

H{l^, 7) = - E "^(^c) Inm(Ae), 

i=l 
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where m(Ac) = / n&AilJ'{x)dx is the invariant measure of Ai. By defining n-th 
refining 7(71) of 7 as: 

n-l 

7'^= U(*iv(x,«))-('=)(7), 

fc=0 

then entropy per unit step of refining is defined by : 

h{iJi, $jv(a;,a),7) = hmn^oo(-^^(Ai, 7)) 

Now, if the size of individual laps of 7(iV) tends to zero as n increases, then 
the above entropy is known as KS-entropy, that is: 

h{n, $iv(a;, a)) = h{ii, ^n{x-, a), 7)- 

KS-entropy, which is a quantitative measure of the rate of information loss 
with the refining, may also be written as: 

r d 

h{iJ,,^N{x,a)) = / iJ,{x)dx\n \ —^N{x,a) |, (V-1) 
J dx 

which is also a statistical mechanical expression for the Lyapunov character- 
istic exponent, that is, the mean divergence rate of two nearby orbits. The 
measurable dynamical system //) is chaotic for /i > and predictive for 
h — As It is proved in Appendix C, the ^n{x, a)) has the following 
form 



h{iJi,^%\x,a)) = In 



^ A^(l + /3 + 2^/;g)^-^ ^ 

■jvi r -'v-i i 



(V-2) 



\ (Sfcio C2kf^^) (Sfc=0 C'2fc+1/^^) / 

Similarly the KS-entropy of coupled chaotic maps can be written as 

roo roo _ d{X\ Ay.) 

^iVi.iVa) = / dXi dX2lJi{x2,X2)hl\ — | . (V-3) 

Jo Jo 0[Xi,X2) 

In Appendix D we have calculated KS-entropy of coupled chaotic maps ^Ni,N2 
for A^i = 2,3,4 and 5, where below we give the results. 
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/ 



+ ln 



iV2(l + /32 + 2v^) 



h{l^, $3,iV2) = In 



+ ln 



2 ' 



\ 



iV2(i + /32 + 2v;^)^^-^ 

-iV2, f iV2 — 1 1 

\ I2^fc=0 '-'2k P2 ) \2^k=0 '-^2fe+lP2 ) I 



X 



$4,Ar2) = In 



21 



\ 



iV2(l + /32 + 2v/^)^^-i 

rJV2l f Af2 ~ 1 1 

, K2^k=0 ^2k P2 ) K2^k=0 ^2k+\P2 ) I 



+ ln 



4(Vl + V^+A/f)'V 



(V-4) 



(V-5) 



X 



21 



(V-6) 



and 
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+ ln 



iV2(l + /52 + 2v^) 



A'2-l 



r 2 1 r 2 ~ 1 

(^2 



5-2V2+ pi, + ^ + 



f32 \ P2 P2 



+ -^V(5-2v^ + /3io)) 



5 + 2 V2 + + ^ + 



2^V^io 2^ 



/?2 N /52 /?2 



+ -f V(5 + 2V2 + /3io)) 



X I (Wl + ^V2 + /5i, + ^ + 



h \ — h — 

/92 \ 02 P2 



+ f^/(l + ^V2 + /30) 



X I (Wl-^V2 + /5io + ^ + 



2^V^io 2^ 
V \ — V — 

P2 \ P2 P2 



ft/(l-^V^ + /?lo)) 



(V-7) 



VI Simulation 

Here in this section we try to calculate Lyapunov characteristic exponent of 
maps {II— 3, 4) and {III— 6, 7) in order to investigate these maps numerically. 
In fact, Lyapunov characteristic exponent is the characteristic exponent of 
the rate of average magnification of the neighborhood of an arbitrary point 
^0 = (^^lo) 2:20) and it is denoted by A(ro) which is written asfl^: 

k k 



./^x -, , , O $ O .... O ($ O $ O .... O <|))2) 

A ro = hm In ^ , — 



1- , dXi{xk,a2) dX2{xk,ai) 
A(ro) = hm 2^ In . (VI-1 
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where = $ o $ o .... o $ . It is obvious that A(ro) < for an attractor, 
A(ro) > for a repeller and A(ro) = for marginal situation. Also the 
Liapunov number is independent of initial point tq, provided that the motion 
inside the invariant manifold is ergodic, thus A(ro) characterizes the invariant 
manifold of $ as a whole. For the values of parameter a, such that the map 



$ be measurable, Birkohf ergodic |15]theorem implies the equality of KS- 



entropy and Liapunov number, that is: 

/i(/^,$) = A(fo,$), (VI-2) 

Comparison of analytically calculated KS-entropy of maps (see Figures 8-11) 
with the corresponding Lyapunov characteristic exponent obtained by the 
simulation, indicates that in chaotic region, these maps are ergodic as Birkohf 
ergodic theorem predicts. In non chaotic region of the parameter, Lyapunov 
characteristic exponent is negative definite, since in this region we have only 
single period fixed points without bifurcation. Therefore, combining the 
analytic discussion of section V with the numerical simulation we deduce 
that, these maps are ergodic in certain values of their parameters as explained 
above and in complementary interval of parameters they have only a single 
period one attractive fixed point, such that in contrary to the most of usual 
one-dimensional one-parameter family of maps they have only bifurcation 
from a period one attractive fixed point to chaotic state or vise versa. A 
comparison of these results with those of obtained in section V, verified the 
accuracy of implemented computation methods. 
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We have given hierarchy of exactly solvable one-parameter families of one- 
dimensional chaotic maps with an invariant measure, that is measurable dy- 
namical system with an interesting property of being either chaotic (proper to 
say ergodic ) or having stable period one fixed point and they bifurcate from 
a stable single periodic state to chaotic one and vice-versa without having 
usual period doubling or period- n-tupling scenario. Also by coupling these 
maps in an appropriate method, we have obtained some coupled chaotic maps 
with invariant measure, where they also bifurcate without period-n-tupling 
scenario to chaos. Since the invariant measure plays an important roles in 
analytic calculation of important equations such as KS-entropy, therefore it 
would be interesting to find the other measurable one parameter maps, many 
parameters, coupled and higher dimensional maps with invariant measure, 
which are under investigation. 



Coupled maps with Invariant Measure 18 

Appendix A: Derivation of measure of one-parameter families of 
chaotic maps 

In order to prove that measure {IV — 3) satisfies equation {IV — 1), it is 
rather convenient to consider the conjugate map. 

l>Ar(a;, a) = ^ tan^(A^arctan a/x) (A-1) 

with measure p,^^ related to the measure /i^^^ with the following relation: 

Denoting ^n{x, ck) on the left hand side of {A — 1) by y and inverting it, we 
get : 

1 I JvTT 

ajfc = tan^(— arctan y'7/Q;2 + — ) k = l,..,N. (A-2) 
Then, taking derivative of Xk with respect to y, we obtain: 

Substituting the above result in equation {IV — 1), we get: 

k 

Now, by considering the following anatz for the invariant measure p,^^ {y) : 

the above equation reduces to: 

which can be written as: 

1 + a'^y a f 3 — l\ d ,^ , ,,, . 
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To evaluate the second term in the right hand side of above formulas we can 
write the equation in the following form: 

— a^y cos^{N arctan -s/x) — sin^(A^ arctan ^/x) 
M ( 



-1 



;i + a;) 



N 



k=0 k=0 



constant 



TV 



+ 

where Xk are the roots of equation {A — 1) and they are given by formula 
[A — 2). Therefore, we have: 

d 



d(3- 



d 



In 



;i - /?-^)^(a27/cos^ (AT arctan^-/?-!) - sin2(iV arctan 



N(3 PN{l + a''y)Aa[ 



+ 



/3-1 ■ 

In deriving of above formulas we have used the following identities 

A{-x) 



(A-7) 



cos(A^ arctan ^/x) 



N_ 1 
2 



sin(A'" arctan ^/x) — \fx 



(1 + x) 

B{-x) 



2 



(A-8) 



inserting the results (A — 7) in (A — 6), we get: 

1 + o?y 1 + o?y 



l + (3y p(^) 

Hence to get the final result we have to choose the parameter a as: 
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Appendix B: Derivation of measure of coupled chaotic maps. 
In order to prove that measure {IV — 7) satisfies equation {IV — 5), it is 
rather convenient to consider the conjugate map. 

Now, denoting Xi by yi and X2 by 1/2 and inverting the conjugate map ^Ni,N2 
we get: 



Xkik2 = tan^ arctan ^fy^a^J^ + ^) 



(B-1) 



= tan^ (-^ arctan \/y2ai + 



2 _|_ fc^ZL^ 
N2 ' 



then taking the derivative of {B — 1) with respect to yi, y^'- 



^3^fc2 Q;2 (l"t"'^fe2) \/^k2 



(B-2) 



9j/2 A^2 ■ (1+053/2) ■ x/y2 

and substituting the above results in equation {IV — 1), we obtain: 



i^{yi,y2) = 

ki,k2 



dxk2 dxkjk2 
dy2 ' dyi 



l-^{-^kik2 ) •^k2 ) 



_ 9l{Xk2) v/3:fcifc2(l + Xk-^kj) tt2 (1 +Xk^) ^fXk2 -., . 



ki,k2 



(B-3) 



Now, considering the following anzatz for invariant measure 

we come to the following equation: 

7/2) = - X] ^^''^ 1 7(2^^2) V- 9'i(^fc2) (l+a^fciifc2) 



(B-4) 



X 



TT^ 9y2 'Vyr'(l + ^i(a;fc2)yi) ^1 ' + Pi{^k2)xkik2y 

(B-5) 



using the relation {III — 3) in the last part of above relation, it reads 
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Nl ' {1 + Pl{Xk2)Xk,k2) BN,{Pl\Xk2)) 



X- 



therefore, the relation {B — 5) reduces to: 

, „ N 1 V- "yj^k^) gi{xk2)ANAPr\xk2)) 
^^'"'^^ --iKV^- B.M\xk2)) 

1 dxk.. 

On the otherhand, we know that the measure fl{yi,y2) is 



or it can be written as 



7(1/2) _ ^ dxk^ -fixk^) 



i(z/2) k2 Mm 



(B-6) 



Kvu 2/2) = -•TT-TSr^fr^' (B-7) 
which is possible if [3i{xk2) and l3i{y2) to be relate with 

and 7(^/2) can be given in terms of ^{xk) as: 

I . ^ ( \/Pi{y2) dxkA , . 



Therefore, T*-^^-* satisfies the PF equation of map $iv, hence it is propor- 
tional to its invariant measure, that is 

7(^2) 1 



,1^^) 7r-(l + /?2t/2)' 



(B-9) 
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Appendix C: Derivation of entropy of one-parameter families of chaotic 
maps 

In order to calculate the KS-entropy of the maps ^n{x, cc), it is rather con- 
venient to consider their conjugate maps given by (// — 2), since it can be 
shown that KS-entropy is a kind of topological invariant, that is, it is pre- 
served under conjugacy map. Hence we have: 

h{ii, = !>). 

The integral — 1) can be written as 

h(iJ,,^N(x,a)) = / fl(x)dxlii. \ — (— . tan^(A^arctan v^)) I 

ax a'' 







we have 

f, N 1 sin Ar(arctan ^/x) A 

using the relation (A — 8), we get 

hi,,^.ix,a))^-l ^^^_^^^3 Ij. (C-1) 

We see that polynomials appearing in the numerator ( denominator ) of 
integrand appearing on the right hand side of equation (C — 1), have ^^^^ 
(■^) simple roots, denoted hy x^ k = 1, [^^-^] {x^ k = 1, [y]) in 
the interval [0, 00) . Hence, we can write the above formula in the following 
form: 

r JV-li 



w ^ / r ^dx , /AT (l + a;)^"inL-i 

h{ii, ^n{x, a)) = - / ' In —A ^ 

71 Jo y'x^l + fJx) {a-^ t-tItJ \ ^ _ \3 



rifcii \x-x^ 

Now making the following change of variable x — ^tan^ ®, and taking into 

account that degree of numerators and denominator are equal for both even 
and odd values of N, we get 
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h(ii, <^>n(x, a)) = - d^{ln(— ) + (iV - 1) In I /5 + 1 + (/3 - 1) cos ^ I 
TT ^0 a"' 

Af-l j I" JVj 

J In I l-xf/3 + (l + xf/3)cos^ I -3^ In I 1 + (1 + cos^ |}. 



k=l 



k=l 



Using the following integrals: 
1 Z-^. 



— / In I a 

TT ^0 



+ 6 cos 9 



In 


a+^a'^ — b'^ 1 
2 1 


1 « > 


b 1 


< 








In 


b 1 
2 1 


1 « l<l 





we get 



h{iJ,, ^n{x, a)) = In 



7V(1 + /3 + 2v;g)^-^ 

\ {J2k=0 C2kP''){J2k=0 ^2k+l(^^) ) 



(C-2) 
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Appendix D: Derivation of entropy of coupled chaotic maps. 
The KS-entropy of ^2,N2 given in {V — 3) can be written as 

roo rco 8X2 

h{fl,^2,N2) = / il{xi)dxi / jl{x2)dx2{ln \ — — | + In | — ^ |). (D-1) 
Jo Jo 0x1 0x2 



Now, using the formulas (/// — 5) the second integral reads 

d , \_ 



foo POO Q \ 

I fl(xi)dxi / jl(x2)dx2ln\ — — (^. tan^(A^2 arctanA/x2)) 
Jo Jo 0x1 ai 



or 



= / iJi[xi)dxi / iJi[x2)dx2hi 

Jo Jo \ a2 v^(l 

where using the relation (74 — 8), we get 



sill .V2 ( f ii'ct an ^^1^) \ 
+ X2) cos^ A^2(arctany^) y' 



^dX2 /iV2 , (1 + X2)'''-'B{-X2) 



xl{A{-X2W 



(D-2) 



IT Jo ^/x^{l + (32X2) V^i 

Finally, calculating the above integral with the prescription of Appendix C, 
we obtain: 

iV2(l+/32 + 2v^)^2-i 



/ 



In 



r N2 1 r JV2 — 1 1 

\ \Z^k=0 ^2k P2 ) {2^k=0 '-^2fc+lP2 ) I 



:d-3) 



Now, using the relation (// — 3) together with the invariance of measure, the 
first integral in (D — 1) can be 



TT^ Jo 



\i{x\)dx\ / ii{x2)dx2\^ 



d P{X2) AliwT)) 



dx, (3{X2y Bl{^) 



X 



sin^(2 arctan ^/xi 



— / jj,(xi)dxi / fl(x2)dx2ln 

TT^ JO Jo 



cos^(2 arctan ^/x^) 

Pix2) (1 + 



f3{X2y A2{P{x2))B2{f3{x2)) 



(D-4) 



Due to the invivariance of measure the integral of the expression -stt^^ van- 



P{X2) 
P(X2) 



ishes, since according to Appendix E, we have 



dxn{x) In {p{X2)) 



dxn{x) In {j3{<^{x))) 



dxfi^x) In [I3{x2)) 
(D-5) 



Coupled maps with Invariant Measure 



25 



Making a change of variable x = ^ tan^ | in the first integral and omitting 
the expression it reduces to 

= - de^iin ((/?2(i + Ao) + g) + (/?2(i + Ao) + g) cose)) 

n Jo 

In {^2 + /92Cose) + In + cos^ + ^ cos2e) + In {A + Bcos^ + Ccos2^)^ 

' ^= 1(3(1 + /3i„) + ^ + 3(fn 
i? = Kl + Ao-(i)') 

, C' = |((l + /30 + ^ + (in, 
The above expression can be calculated by using the following integral: 



Which 



1 



TT Jo 



d9l-a{A + Bcos9 + Ccos29) =21nA 



A 



TT 



1 NA-3C + ^+C)--B- ^ V A + B + C -V A-B + C 



\ 



Y 

/ 

(D-6) 



where the above integral has been evaluated by using the well known mean 
values theorem of analytic function 

- / de\n\f{zo + Rexpi9) \ =| f{zo) \ 

TT Jo 

by choosing f(z) — a + PexpiO + ^exp2i9 . Therefore, the first integral in 
(D-1) reads 
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In- 



Now, combing it with (D — 3) we get the expression ( V-4) . 



2 ' 



(D-7) 
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Appendix E: Proof the integral identity {D - 5) 
Denoting ^{x) by y and inverting it we get 

Xk^^~\y), k^l,2,...,n 

with n as the number of branch, we can write the lefthand side of relation 
{D - 5) as 

/ dxn{x)f{(f){x)) = J2 / dxkPL{xk)fmxk)), (E-1) 

where {xi)k{{xf)k denote the end coordinate of the k-th branch) denote it 
initial coordinate of the k-th branch now, making the change of variable for 

Xk toy using the relations, oo = ^{{xf)k), = ^((a;,)^) and f{y) = f{^{xk)), 
together with the FP equation we obtain 



roo 

= / dy^y)fx{y). 
Jo 



(E-2) 
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Figures Captions 
Fig.l. Plot of $2(2;, a) map for a — 0.19 (dashed curve) and plot of $3(2;, a) 
map for a — 0.19, (solid curve). The location of maxima and minima and 
their values(0 or 1) are independent of parameter a, as it is shown in the 
figure. 

Fig.4. Plot of $4(x, a) map for a — 0.19 (solid curve) and plot of $5(^,0;) 
map for a — 0.19, (dashed curve). The location of maxima and minima and 
their values(0 or 1) are independent of parameter a, as it is shown in the 
figure. 

Fig. 3. Bifurcation diagram of $2(2^,0;), for < a < 2, it is ergodic and for 
< a < cx), it has stable fixed point at x = I.. 

Fig .4. Bifurcation diagram of $3(0;, a), where for | < o; < 3, it is ergodic 
and for < q; < |, it has stable fixed point at x — 0, while for 3 < a < 00, 
it has stable fixed point x — 1. 

Fig. 5. Bifurcation diagram of $2,3, where for < ai < 2 and | < 0:2 < 3 it 
is ergodic, while for 2 < ai < 00 and < q;2 < | (3 < a;2 < 00) it has fixed 
point at [0,0] ([0,1]). 

Fig. 6. Bifurcation diagram of $3.2, where for | < ai < 3 and < q;2 < 2 it is 
ergodic , while for < Q!i < | (3 < ai < 00) and < q;2 < 2 (2 < q;2 < 00) 
it has fixed point at [0, 0] ([1, 0]). 

Fig. 7. Bifurcation diagram of $3,3, where for | < q;i < 3 and | < q;2 < 3 it 
is ergodic, while for < ai < | (3 < ai < cxo ) and < q;2 < | (3 < ^2 < 00 
)it has fixed point at [0, 0] ([1, 1] ). 

Fig. 8. Sohd curve shows the variation of KS-entropy of $2(2;, «), in terms of 
the parameter a dots shows the variation of Lyapunov characteristic expo- 
nent of $2(2;, a), in terms of the parameter a. 
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Fig. 9. Solid curve shows the variation of KS-cntropy of $3(0;, a) in terms of 
the parameter a, dots shows the variation of Lyapunov characteristic expo- 
nent of ^3{x, a), in terms of the parameter a. 

Fig. 10. Solid surface shows the variation of KS-entropy of $2,2 in terms of the 
parameters ai and a2, while dotted surface shows the variation of Lyapunov 
characteristic exponent of $2,2 in terms of the parameters ai and a2- 
Fig. 11. Solid surface shows the variation of KS-entropy of $4^4 in terms of the 
parameters ai and a2, while dotted surface shows the variation of Lyapunov 
characteristic exponent of $4,4 in terms of the parameters ai and a2- 
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